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The reversible control of a single spin of an atom or a molecule is of great interest in Kondo physics and a 
potential application in spin based electronics. Here we demonstrate that the Kondo resonance of 
manganese phthalocyanine molecules on a Au(l 1 1) substrate have been reversibly switched off and on via a 
robust route through attachment and detachment of single hydrogen atom to the magnetic core of the 
molecule. As further revealed by density functional theory calculations, even though the total number of 
electrons of the Mn ion remains almost the same in the process, gaining one single hydrogen atom leads to 
redistribution of charges within 3dorbitals with a reduction of the molecular spin state from S = 3/2 to S = 1 
that directly contributes to the Kondo resonance disappearance. This process is reversed by a local voltage 
pulse or thermal annealing to desorb the hydrogen atom. 

Control over charge and spin states at the single molecule level is crucial not only for a fundamental 
understanding of charge and spin interactions but also represents a prerequisite for development of 
molecular electronics and spintronics 13 . While charge manipulation has been demonstrated by gas 
adsorption 4 and atomic manipulation 5 , the reversible control of a single spin of an atom or a molecule has been 
challenging. Typically, atomic or molecular spin can be probed as a Kondo effect that manifests itself as a 
conductance anomaly at the Fermi level when it is coupled to a metallic system 6 ' 7 . However, an effective method 
to manipulate molecular spin both individually and ensemble has been limited 8 . In this paper, we present the 
reversible control of a single spin in manganese phthalocyanine molecules on a Au(lll) substrate. This process 
can be reversed by a local voltage pulse or thermal annealing to desorb the hydrogen atom, accompanied by a 
recovery of the molecular Kondo resonance. 

Molecules with a single metal ion, including metal phthalocyanines (MPcs), metal porphyrins (MPs) and their 
derivatives, are ideal prototypes to study spin related phenomena due to their versatility and the tunability of 
charge and spin properties 9,10 . When magnetic MPcs or MPs are placed on metal surfaces, the Kondo resonance 
induced by the interaction between a localized spin of the molecule and conduction electrons from the substrate 
manifests itself as a narrow, pronounced peak in the density of states close to the Fermi level. While many attempts 
have been applied to govern the Kondo process, including conformation variation 11 " 15 , alteration of adsorption 
site 16,17 , quantum size effects 18 , ligand attachment 19 , molecular assembly 20 and atomic doping 21 , most are either 
irreversible or cannot manipulate the Kondo effect at both single- molecule scale and ensemble; an effective method 
of molecular spin manipulation has been limited. We show that spin state of a single manganese phthalocyanine 
(MnPc) molecule on the Au(lll) can be reversibly switched by chemical absorption and desorption of a single 
hydrogen atom. Our novel reversible spin control scheme can be easily realized at both ensemble of molecules and 
single- molecule level, which opens up new avenue of broader applications based on molecular spin state. 

Results 

Figure la shows a typical low temperature scanning tunneling microscopy (LT-STM) image of MnPc molecules 
on an Au(l 1 1) substrate, in which every MnPc molecule possesses three inherent unpaired electrons with a total 
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Figure 1 | Topography and dl/dV curves of the MnPc molecules on the Au(lll) before and after hydrogen atom adsorption, (a) Topography of the 
MnPc molecules on the Au( 111) before hydrogen addition, showing molecules with bright center. Scanning bias: -0.2 V; Tunneling current: 10 pA; Scale 
bar, 5 nm. (Inset) High-resolution image of a single molecule. Scale bar, 0.7 nm. (b) Topography of the MnPc molecules after hydrogen atom 
decoration, showing some molecules with suppressed center. Scanning bias: —0.2 V; Tunneling current: 10 pA; Scale bar, 5 nm. (Inset) High-resolution 
image of a single molecule with suppressed center. Scale bar, 0.7 nm. (c) Sequential variation of the dl/dV spectra recorded at the center of a MnPc 
molecule induced by the absorption and desorption of a single hydrogen atom. For comparison purpose spectra were vertically shifted 0.1 nA/V and 
0.33 nA/V for middle and upper curves, respectively, (d) dl/dV spectra acquired at the center of the MnPc molecule under external magnetic field at 0.4 K, 
showing the splitting of Kondo resonance. For clarity the successive spectra were vertically shifted by 0. 1 nA/V. 



spin of S = 3/2. After deposition of 0.05 monolayer of MnPc onto a 
22 X ^3 herringbone reconstructed Au(lll) surface, isolated MnPc 
molecules preferentially occupy the elbow and face centered cubic 
(fee) sites. The STM image unambiguously confirms an intact 
molecular structure after thermal deposition with exhibition of a 
protruding "cross" feature with four-fold symmetry 1618 . After the 
MnPc/Au(lll) sample was dosed with 30 Langmuir (L) of H 2 at 
room temperature (RT) and subsequently cooled down to 4.2 K, 
all molecules retained the cross feature but some showed a depres- 
sion at the center. After extensive H 2 dosage, all of the original MnPc 
molecules can be turned into molecules with a depression at the 
center (Fig. lb and Supplementary Fig. 1). 

The variation of the molecular morphology after the H 2 gas dosage 
can be attributed to the chemical adsorption of a single hydrogen 
atom onto the MnPc molecule, according to key control experiments 
summarized as follows: (1) This change of molecular morphology is 
reversible. We observed that by applying a positive voltage pulse (e.g. 
2.0 V for one second with feedback loop off) the bright protrusion 
feature of the MnPc molecules can be recovered on a depression. This 
behavior is consistent with a tip induced detachment 4 of adsorbates 
from the MnPc molecule. Alternatively, the depression feature can be 
switched back to a bright protrusion by thermal annealing of samples 
at 500 K for 10 mins. (2) The apparent height of the molecule re- 
mains unchanged when the H 2 is dosed at 4.2 K instead of RT, even 
with very high dosage capacity (1000 L). We observed that intact 
H 2 molecules can only reside around the MnPc molecules without 
modification of the intrinsic molecular structure (Supplementary 
Fig. 2). This suggests that the interaction between the H 2 and the 
MnPc molecules is weak and absorption of molecular hydrogen 
cannot be the origin of the process demonstrated in the Fig. 1. 
(3) A threshold voltage of +1.28 V is required in order to detach 
the adsorbate from the host MnPc molecule (Supplementary Fig. 3). 
This threshold voltage suggests the binding between adsorbate and 
the molecule is chemical adsorption rather than physical adsorption, 
which demonstrates that the absorption of the H 2 cannot be the 
origin of the process shown in the Fig. 1. (4) Lastly we also carried 
out one control experiment by direct dosage of atomic hydrogen 
gas 22 (Supplementary Section 2.), and similar features of the 



molecular morphology change were also obtained. From these key 
experiments we conclude that the depression feature after dosage 
of hydrogen gas is due to adsorption of atomic hydrogen (labeled 
as "H-MnPc"). 

We explored and compared in detail the difference of electronic 
structure between the MnPc and the H-MnPc states by measuring 
differential conductance (dl/dV) spectra. As shown in the upper 
red curve in Fig. lc, the dl/dV spectra acquired at the center of the 
original MnPc molecules consistently show a pronounced step 
shaped feature at zero bias, regardless of the molecular adsorption 
configurations and locations on the Au terraces. Such step shaped 
feature in dl/dV spectra is very sharp and close to the Fermi level. 
This spectra feature manifests linear splitting at the presence of mag- 
netic field (Fig. Id) as well as temperature dependence (Supplemen- 
tary Fig. 8), and thus it can be safely attributed to the Kondo 
resonance 6 ' 7 ' 23 ' 24 . In contrast, the dl/dV spectra of the H-MnPc state 
(middle blue curve in Fig. lc) are featureless in this energy range. 
Importantly, once the H-MnPc state was switched back to the MnPc 
state, its electronic structure including the Kondo resonance can be 
fully recovered (bottom red curve in Fig. lc) in addition to the recov- 
ery of the topographic feature. 

Further insight into the reversible switching of the Kondo res- 
onance between the MnPc and the H-MnPc states can be gained 
by density functional theory (DFT) calculations (see Method for 
more details). Topographic features of both the MnPc and H- 
MnPc states can be well reproduced by STM simulations by integ- 
rating electron density from —0.4 eV to the Fermi level according to 
the Tersoff-Hamann approximation 25 . The projected density of 
states (PDOS) shows that the bright protrusion at the molecular 
center of MnPc is due to the d-orbital character of the Mn ion near 
the Fermi level (Fig. 2a, and Supplementary Fig. 5). The calculated 
binding energy of 1.75 eV between the hydrogen atom and MnPc/ 
Au(lll) indicates that H on MnPc/Au(lll) is of a chemical adsorp- 
tion with binding energy in the same range of our experimental 
value. Therefore, both experimental results and theoretical calcula- 
tions suggest that attaching and detaching a single hydrogen atom 
from the central Mn ion is the physical origin of spin switching 
process. 
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Figure 2 | STM simulation and PDOS of the MnPc/Au(lll) before and 
after hydrogen adsorption (a) STM simulation of the MnPc/Au(lll) 
and the H-MnPc/Au(lll), confirming the topographic feature in Fig. la 
and b. (b) and (c) PDOS of d xz > d yz and d z i orbitals of Mn ion in the MnPc/ 
Au( 111) and the H-MnPc/Au( 111), respectively, revealing that spin at the 
Mn ion is reduced from 3/2 to 1. (d) Schematics of reversible control of 
molecular Kondo effect with adsorption (Kondo OFF) and desorption of 
hydrogen atom (Kondo ON). 

Discussions 

An analysis of the electronic structure of MnPc and H-MnPc reveals 
that the number of d- orbital electrons of the Mn ion remains almost 
unchanged. Thus the additional electron contributed by the hydro- 
gen atom is donated to the Pc skeleton. However, the effective 
charges in individual d orbitals are redistributed. Before hydrogen 
adsorption, the Mn ion in the MnPc molecule has a net spin of S = 
3/2. This is due to two partially paired orbitals, d xz and d yz , and two 
unpaired orbitals, d xy and d z i (Fig. 2b). The Kondo resonance 
observed in the dl/dV spectra in this case originates from an 
exchange interaction between the localized spin of the magnetic 
Mn ion and the conducting electrons of the Au(lll) substrate 26 . 
After atomic hydrogen decoration (i.e., in the H-MnPc state), the 
spin polarization decreases in the d z i orbital: Before adsorption the 
spin-up charge in this state is close to one while the spin-down charge 
is close to zero. After adsorption the charge in both spin states is close 
to 1/2. This leads to a decrease of the net spin of the Mn ion from S = 
3/2 to S = 1 (Fig. 2c). Because the Kondo effect strongly depends on 
the coupling between the localized spin and the conducting electrons, 
a reduced molecular spin state directly contributes to the suppression 
of the Kondo effect (Supplementary Section 4) 27,28 . A schematics of 
manipulation of the molecular Kondo effect by hydrogen adsorption 
and desorption is illustrated in Fig. 2d. 

In principle, the appearance or disappearance of the Kondo res- 
onance on this magnetic molecule can be seen as a single bit of 
information. This may have implications for information recording 
and storage at the ultimate molecular limit 10 " 12 ' 29 " 31 , when combined 
with STM capability of atomic manipulation of molecular adsorption 
and desorption. We explore this avenue by investigating an ordered 
molecular array. Fig. 3a shows topography and dl/dV mapping 
acquired simultaneously from a 3 X 4 molecular array. All the 
MnPc molecules within the array are initialized by converting them 
to the H-MnPc state through the adsorption of hydrogen atoms, with 
evidence of featureless dl/dV mapping (i.e. absence of the Kondo 



resonance), while the four-lobe topographic structure of individual 
molecules remains (Fig. 3a). By applying voltage pulses on selected 
H-MnPc molecules, we can precisely address an individual molecule 
within the array and the H-MnPc molecules can be converted back to 
the MnPc state one-by-one. As a result, a pattern of confined dark dot 
feature corresponding to the Kondo resonance of the MnPc mole- 
cules can be created in the dl/dV mapping (Fig. 3b). Importantly, this 
process can be extended to an ultimate close-packed array with the 
highest molecular density (Supplementary Fig. 6 and Supplementary 
Fig. 7), with one example shown in the Fig. 3c. While the inter- 
molecular spacing is as small as 1.4 nm in a close-packed array, all 
the dark-dot features in the dl/dV mapping remain localized at the 
center of the MnPc molecules, and the dl/dV spectra of the Kondo 
switching recorded in the molecular assembly show no difference 
from our observation of isolated molecules (Supplementary Fig. 7), 
regarding characteristic parameters such as the Fano factor (q) and 
the Kondo temperature (T K ). Correspondingly, no change occurs in 
the molecular framework during the spin manipulations within the 
molecular assembly, in contrast to previously investigated systems, 
which involve conformational changes at the periphery of the mole- 
cule 11 ' 1314 . A storage system based on the molecular spin states thus 
describes the ultimate limit of information storage. Furthermore, all 
the MnPc molecules in the arrays exhibit consistent Kondo features, 
suggesting a site-independent Kondo resonance and the inter-par- 
ticle interactions can be weak enough to become negligible in the 
process. (Fig. 1 and Supplementary Fig. 7). Even more importantly, 
such molecular spin switching can be consistently achieved back- 
and-forth for many times with no observable change of the Kondo 
features. To substantiate the robustness of the spin switching process, 
we have quantified both T K and q of MnPc state from a Fano fitting 
(Supplementary Fig. 8) and monitored both parameters during the 
spin switching process by repeatedly absorption and desorption of 
H-atom, with results summarized in the Fig. 3d. That both, the T K 
and q values in the spin manipulation cycles remain constant sup- 
ports the robustness of this process. Atomic localization feature of 
the Kondo resonance switching makes it possible to create more 
complex patterns based on molecular spin state, including one 
example of a digitalization from (000) to (1 1 1) that can be fabricated 
from a larger 3X8 molecular array (Supplementary Fig. 9). 

The ability to reversibly tune the Kondo effect at the molecule level 
by gaining or losing a single hydrogen atom creates a new class of 
simplest molecular Kondo system in which the spin screen can be 
manipulated directly through the control of the local atomic envir- 
onment. Our demonstration of reversible molecular spin control by 
chemical stimuli is different from a few recent works in that it offers 
an in -situ single molecule characterization and control instead of 
ensemble measurement 32 and selection of noble metallic Au sub- 
strate can further avoid complex spin texture of substrate due to 
Rashba type spin-split pair surface states of a semimetal of such as 
Bi 33,34 . Even more importantly, our scheme can allow not only 
ensemble spin manipulation (such as Kondo ON/OFF switching of 
all molecules by thermal annealing and H 2 gas dosage, respectively) 
but also control at the single-molecule level by atomic manipulation 
of single H atom that has not been achieved so far 21 . It thus opens up a 
variety of avenues. Our approach can be readily extended to more 
complex magnetic molecules. The atomically precise engineering 
of molecular spin states provides a unique test ground for probing 
magnetic interactions at the single-molecule level. For example, the 
magnetism in transition -metal ions based magnetic molecules arises 
from unpaired spins residing in the d- orbitals. Controlled coupling 
of its net spins with another unpaired spin of a well-defined hydrogen 
atom provides the simplest conceivable systems potentially exhib- 
iting quantum critical behavior. Our demonstration of robust Kondo 
switching in a molecular array could also be applied as nonvolatile 
classical information recording and storage at the ultimate molecule 
limit without the requirement of molecular conformation change. 
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Figure 3 | Reversible spin switching in molecular arrays, (a) The dl/dV mapping (upper panel) and simultaneously acquired topography 
(lower panel) of a H-MnPc molecular array (Kondo OFF state), (b) The dl/dV mapping (upper panel) and simultaneously acquired topography (lower 
panel) of the same molecular array as a but with a few selected molecules converted to the MnPc state (Kondo ON). For both (a) and (b), the dl/dV 
mappings were taken at 6 mV. Note that the MnPc state (Kondo ON) can be erased and converted back to H-MnPc state (Kondo OFF) by adsorption of 
additional hydrogen, (c) The dl/dV mapping (upper panel) and topography (lower panel) of a close-packed molecular array with pre-designed Kondo 
pattern by switching selected molecules from the H-MnPc to the MnPc states within the array. The <il/<i Vmapping was taken at 6 mV. (d) The Fano factor 
q and Kondo temperature T K of Kondo ON states after multiple- cycles of spin switching, indicating the robustness of the spin manipulation process. 



Future design of more complex molecular framework might offer a 
way to enable and tailor coherent spin coupling through such as 
inter-molecular bridge 35 . 

Methods 

Experimental method. Experiments were carried out in an ultrahigh vacuum (base 
pressure of 1 X 10" 10 mbar) LT-STM system (Unisoku), equipped with standard 
surface processing facilities. An atomically flat Au(lll) surface was prepared by 
repeated cycles of sputtering with argon ions and annealing at 800 K. Commercial 
MnPc molecules (Sigma- Aldrich, 97% purity) were sublimated from a Knudsen-type 
evaporator after thermal purification, while the Au(lll) substrate was held at room 
temperature. One monolayer refers to a close-packed MnPc layer covering the whole 
Au(lll) surface, as examined by STM (Supplementary Fig. 6). STM images were 
acquired in the constant-current mode with bias voltage applied to the sample. 
Spectroscopic measurements were performed by using a lock-in technique with a 
0.5 mV rms sinusoidal modulation at a frequency of 973 Hz. Most images and spectra 
were obtained with electrochemically etched tungsten tips at 4.2 K, but the magnetic 
field dependent measurement was performed at 0.4 K. Measurements of local 
spectroscopy were electronically calibrated by performing dl/dV measurements on 
clean Au(l 1 1) before and after measurements on molecules, based on the well-known 
Au(l 1 1) surface state 36 . The atomic H is generated by cracking H 2 molecules with hot 
tungsten filament, which is placed very close (~5 cm) to the sample. 

Calculation details. The DFT calculations were performed using Vienna ab-initio 
simulation package (VASP) 37 . Perdew- Wang 9 1 exchange correlation functionals and 
the projector augmented-wave method were employed 38 . The energy cutoff for 
plane- wave basis set was 400 eV. A c(7 X 8) supercell containing three layers of gold 



atoms was employed to model an isolated molecule on Au(lll) substrate. To ensure 
that interactions between the periodic slabs through vacuum can be negligible, the 
slabs were separated by a vacuum gap of 20 A. For geometry optimization, the bottom 
layer was fixed, while the adsorbates and the other metal layers were allowed to relax 
until the forces applied to the relaxed atoms were less than 0.01 eV/A. A single T 
point was used in sampling the Brillouin zone due to the numerical limitations (See 
Supplementary Section 8). 
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